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摘 要 :精准 的 预测 评估 是 制定 岩 爆 防 治 策略 的 重要 前 提 。 采 用 FISH 语言 编程 对 FLAC” 数 值 模 
拟 软 件 进 行 二 次 开发 ,基于 动力 数值 分 析 , 将 6 种 岩 爆 预测 评价 指标 引用 到 锦 屏 下 级 水 电站 4 引水 
隧洞 数值 模拟 岩 爆 案 例 中 ,对 岩 爆 的 孕育 过 程 和 形成 机 制 进 行 深入 剖析 。 研 究 结 果 表 明 : 隧 洞 掌 子 
面 右 侧 边 墙 及 拱 腰 部 位 发 生 了 以 剪 切 破坏 为 主 的 突 发 性 破坏 , 洞 室 周围 岩 体 破裂 演化 的 平缓 加 速 
历程 是 围 岩 内 部 能 量 逐 渐 积聚 、 耗 散 的 过 程 ; 在 进行 岩 爆 发 生 位 置 预测 时 , Turchaninov 判 据 较 
Hoek , Russenes 判 据 与 现场 情况 更 为 符合 ,直接 验证 了 其 作为 岩 爆发 生 位 置 判 据 的 合理 性 。 研 究 成 
、 采 可 为 类 似 工程 提供 借鉴 和 参考 。 
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Abstract: Accurate prediction assessment is an important premise for formulating rockburst prevention 
strategies. FLAC?” numerical simulation software was secondarily developed by using FISH language pro- 
gramming, and based on dynamic numerical analysis , after a detailed analysis of the evolution process and 
formation mechanism of rockburst, six kinds of rockburst prediction assessment indexes are applied to No. 
4 diversion tunnel of Jinping II Hydropower Station. The results showed that the shear failure occurred 
mainly in the right side wall and arch haunch of the tunnel face , and the gentle acceleration process of rock 
fracture evolution around the cavern is also the process of energy accumulation and dissipation in the sur- 
rounding rock. When predicting the location of rockburst , Turchaninov criterion is more consistent with the 
field situation than Hoek and Russenes criterion, because it directly verifies its rationality as the location 
criterion of rockburst. The research achievements can provide some reference for similar projects. 


Key words: tunnel ; rockburst ; rockburst prediction; dynamic calculation 


收 稿 日 期 :2020-09-07 修 回 日 期 :2021-10-17 
基金 项 目 :国家 自然 科学 基金 ( No. 51474097;51778215;U1810203 ) ;河南 理工 大 学 青年 骨干 教师 资助 计划 (No. 2017XQC-08) 
通信 作者 : 郭 佳 奇 ,副教授 。E-mail:gjq519@ 163. com 
引用 格式 : 孙 飞 跃 , 刘 希 亮 , 郭 佳 奇 ,等 . 岩 爆 预测 评估 方法 的 动力 数值 分 析 [ 本: 应 用 力学 学 报 ,2022 ,39 (1) :26-34. 

SUN Feiyue, LIU Xiliang, GUO Jiaqi, et al. Dynamic numerical calculation analysis of rockburst prediction assessment methods[ J ]. Chinese 


journal of applied mechanics ,2022 ,39 ( 1) :26-34. 


Á 21 世纪 以 来 ,我 国 在 交通 运输 和 水 利水 电 等 
领域 的 地 下 洞 室 工程 建设 竞相 发 展 ,诸如 西 电 东 送 、 
西 气 东 输 和 南水北调 等 重大 工程 取得 了 重大 进展 。 
由 于 深部 工程 开 挖 诱导 产生 的 岩 爆 动力 灾害 对 地 下 
工程 具有 极 大 的 破坏 性 ,已 造成 施工 设备 损毁 和 施 
工人 员 伤 亡 的 严重 后 果 , 并 产生 了 巨大 经 济 损 
失 "] 。 锦 屏 开 级 水 电站 引水 隧洞 是 我 国 水 利水 电 
工程 中 具有 代表 性 的 案例 之 一 , 最 大 埋 深 达 
2525 m 最 大 地 应 力 高 达 60 MPa ,在 施工 期 间 均 
发 生 了 不 同 程度 的 兰 爆 动力 灾害 。 为 了 给 深 埋 地 下 
工程 的 设计 和 施工 提供 更 符合 工程 实际 条 件 的 科学 
依据 ,和 安全 合理 的 理论 技术 支持 ,对 岩 爆 动力 灾害 
的 预测 评估 方法 进行 系统 地 研究 ,已 成 为 一 种 必然 
地 发 展 趋势 ,并 具有 重要 的 工程 意义 。 

汪 现 如 今 , 岩 爆 动 力 地 质 灾害 的 预测 已 成 为 公认 
的 必 界 性 难题 之 一 571 。 自 岩 爆 出 现 以 来 就 一 直 受 
到 国内 外 岩石 力 学 工作 者 和 工程 研究 人 员 的 高 度 重 
裕 更 是 近 些 年 国内 外 有 关 专 家 与 学 者 共同 关注 的 
HU. Hi. PI Ae CERCA ERROR ,数值 
模拟 和 试验 方面 ,针对 着 爆 动 力 灾害 的 预测 预报 问 
题 焉 展 了 大 量 研究 。 在 预测 预报 评估 理论 研究 方 
EEDI SEINE . 岩 体 失 稳 理论 .突变 理论 、 
能 量 比 法 、 冲 击 能 量 指标 ,应 变 能 储存 指标 和 脆性 指 
数 等 数 十 种 。 在 试验 研究 方面 ,文献 [15] 得 到 
下 省 爆 陆 测 智 能 分 类 模型 。 文 献 [16] 为 准确 预测 
各 防御 时 滞 型 岩 爆 的 发 生 , 对 岩石 失效 的 时 间 行 为 
和 硬 澡 型 岩 爆 的 形成 机 制 进行 了 试验 研究 。 文 南 
[] 四 提出 了 投影 寻 踪 与 董 火 虫 算法 相 结合 的 崖 爆 
预测 评价 方法 。 文 献 [18] 以 室内 试验 为 研究 手段 ， 
Telf AS ERBUR AERE E, E T 2 TRA TOP- 
SIS 法 的 岩 爆 预测 模型 。 在 数值 模拟 研究 方面 ,国内 
外 众多 专家 学 者 开展 了 大 量 的 研究 和 探讨 。 文 献 
[19] 为 了 精准 预测 岩 爆 的 发 生 , 提 出 了 一 种 用 于 模 
拟 KCC 模型 校准 的 方法 。 文 献 [20] 以 试验 和 数值 
模拟 为 研究 手段 ,确定 了 JH-2 模型 参数 ,同时 进行 
了 岩 爆 分 级 预测 。 文 献 [21] 通过 采用 三 维 离散 元 
数值 模拟 软件 ,针对 工程 埋 深 大 、 地 应 力 高 和 地 质 条 
件 复杂 等 问题 ,进行 了 高 地 应 力 条 件 下 超 长 深 埋 隧 
洞 岩 爆 预测 研究 。 

如 前 所 述 ,众多 学 者 从 各 个 角度 对 兰 爆 预报 预 
测 开展 了 深入 而 广泛 的 研究 ,对 其 特征 和 规律 有 了 
一 定 的 认识 和 理解 。 但 预测 评价 方法 过 于 单一 ,未 
能 全 面 考虑 岩 爆 控 制 影响 因素 ,如 :应 力 条 件 与 围 岩 
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更 动态 破坏 之 间 的 关系 , 深 埋 地 下 洞 室 开 挖 过 程 中 
关键 部 位 的 应 力 状态 与 岩 爆 之 间 的 相关 性 。 进 一 步 
了 解 到 在 上 述 预 测评 估 方 法 中 ,因数 值 模拟 法 具有 
经 济 性 、 可 重复 性 以 及 定量 性 等 优势 ,使 得 该 方法 在 
岩 爆 动力 地 质 灾害 的 预测 评估 中 使 用 较 多 。 此 外 ， 
工程 实践 表明 , 围 岩 体内 的 应 力 分 布 状 态 和 应 力 集 
中 现象 对 岩 爆 的 发 生 具 有 严重 的 影响 。 一 般 而 言 ， 
弹性 应 变 能 通常 积聚 在 应 力 较 高 的 区 域 ,其 发 生 的 
岩 爆 可 能 性 也 越 高 。 而 数值 模拟 法 正 是 通过 确定 工 
程 区 域 不 同 部 位 ` 不同 开 挖 阶段 的 应 力 集中 程度 及 
分 布 规律 ,对 任意 地 点 ,尤其 是 未 施工 区 域 ,提前 进 
行 岩 爆 的 危险 性 预测 ,以 及 确定 洞 室 开 挖 中 各 部 位 
出 现 最 大 应 力 的 部 位 和 时 间 , 从 而 对 诱导 岩 爆 的 重 
点 区 域 进行 防治 。 

在 运用 数值 模拟 法 对 岩 爆 进行 预测 时 , 以往 研 
究 大 多 采用 静 力 数值 计算 法 ,然而 此 类 方法 无 法 较 
为 真实 地 反映 岩 爆 的 孕育 演化 过 程 及 其 能 量 动态 响 
应 规律 。 鉴 于 此 ,本 文 对 有 限 差分 软件 FLAC” 进 行 
二 次 开发 ,并 基于 动力 数值 分 析 , 以 实际 工程 案例 为 
依托 ,从 常用 的 6 种 宕 爆 预 测评 价 指标 与 评估 方法 
( Hoek 判 据 、Russenes 判 据 、Turchaninov 判 据 、ERR 
指标 .LERR 指标 和 围 宕 强度 比 ) ,对 岩 爆 的 倾向 性 
进行 预测 评估 ;并 探讨 深部 地 下 洞 室 的 能 量 演 化 规 
律 及 其 动态 响应 ,研究 成 果 可 为 类 似 工程 提供 一 定 
的 参考 依据 。 


1 工程 概况 


锦 屏 工 级 水 电站 4” 引水 隧洞 由 东 往 西 开 挖 , 当 
开 控 至 K9 +728 标 段 时 ,K9 +742 ~ K9 +766 标 段 南 
侧 边 墙 至 拱 脚 部 位 发 生 了 极 强 岩 爆 “( 见 图 1)。 
岩 爆 坑 深度 大 于 2 m, 岩 爆 的 剧烈 冲击 使 锚 杆 拉 断 
且 严 重 扭曲 变形 ”1 ( 见 图 2)。 通 过 现场 勘查 未 发 
现 此 标 段 有 控制 性 结构 面 , 且 围 岩 新 鲜 完整 , 该 洞 段 
围 岩 体 主要 为 T ,大 理 岩 ,4 引水 隧洞 断面 尺寸 见 
图 3 7^ ,依据 现场 监测 的 地 应 力 反 演 结果 '”! ,该 洞 
段 地 应 力 水 平 较 高 ,具体 地 应 力 状态 见 表 1。 

表 1 4 引水 隧洞 岩 爆 段 地 应 力 状态 


Tab.1 The in-situ stress state of 4^ headrace tunnel 


埋 深 / ol oy/ [A TA Tu Ta^ 


m MPa MPa MPa MPa MPa MPa 


1900 -49.81 -51.68 -58.09 -15 =1,23 7:17 
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W 3# 引 水 隧洞 


pum 所 过 方向 E 


4# 引 水 隧洞 K9+728 K9+801 K9+810 


E 修 支 洞 


南 侧 边 墙 至 拱 脚 处 极 强 岩 爆 «mec EESDRIHZ N/A) L3 03H20 


图 1 4 引水 隧洞 岩 爆发 生 位 置 示意 


Fig.1 Rockburst location of 4* headrace tunnel 


图 3 4 引水 隧洞 断面 尺寸 (单位 :m) 


Fig.3 Dimension of 4* headrace tunnel ( unit;m) 


2 岩 爆 预测 理论 及 其 指标 


目前 各 国学 者 在 对 动力 灾害 形成 机 理 进 行 研 究 
时 ,从 多 角度 、 全 方位 、 多 形式 相继 提出 了 一 系列 的 
经 典 理论 ,如 :刚度 理论 .强度 理论 .能 量 理论 、 变 形 
失 稳 理论 和 冲击 倾向 理论 等 “1 ,并 基于 各 自 假设 
提出 了 相应 的 预测 评价 指标 。 在 数值 计算 中 ,通过 
对 三 维 有 限 差分 软件 进行 二 次 开发 ,采用 FISH 语言 
编制 数据 接口 程序 来 实现 式 (1) ~ 式 (6) 的 表达 ,并 
监测 所 有 计算 单元 的 变化 过 程 。 


2.1 Hoek 判 据 


O mx/ 0c =0.34( 少 量 片 帮 , | 级) 

m. 0c =0.42( 严 重 片 帮 , IL SR) 

m, Oc =0.56( 需 重型 文 护 , 亚 级 ) 

O na/Cc >0.70( 严 重 岩 爆 ,区 级 ) 
式 中 :ou 为 隧洞 断面 最 大 切 向 应 力 ;ore 为 岩石 单 轴 
抗 压 强度 。 


(1) 


2.2 Russenes 判 据 


a ac <0.2( 无 岩 爆 活动 ) 
op/oc 20.2 ~0.3( 低 等 岩 爆 ) 
os/oc =0.3 ~0.55( 中 等 岩 爆 ) 
oo/oc >0.55( 严重 岩 爆 ) 
式 中 :0 为 洞 室 切 向 应 力 ;oc 为 岩石 单 轴 抗 压强 度 。 


2.3 Turchaninov 判 据 


(2) 


c, ta, /a. <0.3( 无 岩 爆 活动 ) 
gg+oi/oc =0.3 ~0.5( 有 岩 爆 可 能 (3) 
0,*30,/0, 20.5 -0.8( —EZ REG) 
c,*o,/0, >0.8( 有 严重 告 爆 ) 
式 中 :0 为 洞 室 切 向 应 力 ;or 为 洞 室 轴 向 应 力 ;oc 
为 岩石 单 轴 抗 压强 度 。 


2.4 ”能 量 释放 率 ( ERR) 


W, = L(s Edi 


m^i^i 


| (4) 
ERR = W/V 


sp, 为 某 一 动态 开 挖 步骤 中 国 岩 所 释放 的 能 
1858, 为 本 步 开 挖 而 暴露 出 的 表面 积 ;u; 为 开 挖 引 
起 的 围 宕 次 生 位 移 ;7; 为 本 步 开 挖 前 围 岩 中 的 表面 
牵引 力 ;T 为 开 挖 岩 体 总 体积 。 


2.5 局 部 能 量 释放 率 (LERR ) 


LERR = Uma — U, 

U imax = lo +0, +0; -W( 0,0, t 0,0, +0103)] 

U, is = lo? +0% +0% -2v(0',0, +0',0', +0'10')] 
(5) 


式 中 :Ui 为 第 i 个 单元 脆性 破坏 前 的 弹性 应 变 能 密 
度 最 大 值 ; Ui, 为 第 i 个 单元 脆性 破坏 前 的 弹性 应 变 


能 密度 最 小 值 ;oi os .os 为 单元 应 变 能 最 大 值 对 应 
的 3 个 主 应 力 ;c oa 为 单元 应 变 能 最 小 值 对 
应 的 3 个 主 应 力 ;5 为 弹性 模 量 ;v 为 泊 松 比 。 

2.6 围 岩 强度 比 


max 《<4( 可 能 发 生 严 重 岩 爆 ) 

max 74 ~7( 可 能 发 生 轻 微 或 中 等 岩 爆 ) (6) 
maxs >7( 无 岩 爆 ) 

式 中 :oc 为 岩石 单 轴 抗 压强 度 ; 0 为 图 岩 最 大 地 
应 力 。 


3 ” 岩 爆 预测 的 数值 模型 


[on 
Ce 


Oo/O 


3L 数值 模型 及 边界 约束 条 件 


为 保证 数值 模拟 结果 真实 可 靠 ,在 采用 三 维 有 


限 次 分 软件 进行 求解 动力 学 问题 时 ,依据 圣 维 南 原 
理 和 隧洞 开 挖 的 影响 范围 局 。 建 立 的 数值 模型 

为 200 m,Z Jr [8] Jj 100 m, Y 7r [3] 2j 50 m, 数 值 
模型 共 1 050 500 个 六 面体 单元 1 113 381 个 网 格 节 
Je ROS TL ES] 4 ,监测 点 的 布置 见 图 5。 在 

计算 时 ,阻尼 设置 为 Rayleigh 阻尼 、 最 小 临界 阻 
JEMI 0. 04 .最 小 中 心 频率 取 550 Hz。 该 数值 模型 Z 
方 鸟 的 上 边界 为 应 力 边界 条 件 ,施加 58.09 MPa. 垂直 
Apt Z 方向 的 下 边界 了 方向 的 前 后 边界 和 区 方向 
的 轼 右边 界 受 位 移 约束 ;为 了 减 小 或 消除 模拟 计算 
所 窟 生 的 弹性 波 反 射 , 在 模型 外 围 边界 设 定 为 静态 
边疆 ,以 期 为 计算 模型 提供 等 同 于 无 限 场地 的 约束 
AR ,如 图 6 所 示 。 


4 引水 隧洞 


排水 洞 


4 
Yox 


图 4 数值 模型 
Fig.4 Numerical model 


3.2. ”爆破 荷载 作用 形式 


在 数值 分 析 中 爆破 丛 载 常 假定 为 一 个 三 角形 冲 
击 波 ”” ,如 图 7 所 示 。 三 角形 函数 爆炸 荷载 历程 
曲线 表达 式 , 如 式 (7) 。 通 过 对 FLAC” 软 件 进行 二 
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次 开发 ,采用 FISH 编写 程序 语言 来 进行 爆破 荷载 的 
施加 (去 为 0.3 ms, t, 29 1.5 ms, KI fag EIE (EE N 


60 MPa ,并 利用 APPLY 命令 施加 到 隧道 开 挖 轮 
廓 面 ) o 


图 5 4 引水 隧洞 监测 点 位 置 


Fig.5 Monitoring point position of 4* headrace tunnel 


图 6 自由 场 边界 示意 
Fig.6 Free field boundary 


图 7 爆破 荷载 曲线 
Fig.7  Blasting load curve 


0 t«0,t»t, 


=P; O0 xt «t, 


P(t) =; t (7) 


IP IPOD 为 任 一 时 刻 的 爆破 荷载 压力 值 ; Po 为 爆 
破 荷 载 峰值 ;i, 为 爆破 荷载 上 升 至 峰值 的 时 间 ;i 为 
爆炸 荷载 的 正 压 作用 时 间 。 


3.3 ”本 构 模 型 及 材料 力学 参数 
在 数值 计算 中 ,数值 模型 的 本 构 关系 采用 描述 
硬 岩 力学 行为 的 弹 塑 性 模型 “1 。 围 兰 的 物理 力学 


参数 参照 文献 [22 ] ,如 表 2 所 示 。 本 模拟 对 岩石 岩 
性 作出 假设 :岩石 为 均 质 各 向 同性 的 连续 体 ,符合 


30 
弹 塑性 强度 准则 , 材料 参数 满足 弹 塑性 本 构 模 型 
关系 |。 
表 2 岩 体 的 物理 力学 参数 


Tab.2 Physical and mechanical parameters of rock 


ME yp WRD WEU BRN BRA 
模 量 / 有 峰值 / ”残余 值 /” 初始 值 / 峰值 / 
CPa MPa MPa (9) (95 


27.62 0.256 34.36 9.87 29.93 39.93 29.20 


3.4 计算 分 析 与 岩 爆 预 测 


由 图 1 可知,44 引 水 隧洞 在 洞 室 掌 子 面 后 方 约 
80 m 范围 内 发 生 了 岩 爆 .| 。 该 隧洞 开 挖 为 上 下 人 台 
阶 法 , 仅 对 洞 室 进行 上 台阶 开 挖 时 ( 开 挖 高 度 8.5 m, 
见 图 3) ,4 引水 障 洞 就 发 和 了 岩 爆 现象 。 因 此 ,在 
对 该 隧洞 进行 岩 爆 倾向 性 分 析 时 ,选取 了 发 生 岩 曝 
的 中 间 位 置 ( 标 段 为 K9 +765 附近 ) 。 
ONIL) 基于 Hoek 判 据 分 析 
通过 数值 计算 获得 4* 引 水 隧洞 横 剖 面 的 Hoek 
ee 8 .图 9 
D 


o 


oO oOoooooooo 
—i3U 4 UO 100 


图 8 Hoek 岩 爆 判 据 分 布 
Fig.8 Hoek rockburst criterion distribution 

由 图 8 可 知 ,采用 Hoek 判 据 分 析 时 ,隧洞 掌 子 
面 的 左 侧 边 墙 与 左 拱 脚 附近 处 最 大 指标 值 约 为 
0. 9 ,指标 值 大 于 0.7 的 深度 约 为 1.5 m; 右 侧 边 墙 与 
右 拱 脚 附近 处 最 大 指标 值 约 为 0. 85 ,指标 值 大 于 
0.7 的 深度 约 为 0.8 m, Bil] TEJT ZUG , BED ZU 
PER 、 边 墙 及 拱 脚 处 出 现 弹 性 应 变 能 积聚 现象 ,最 大 
积聚 深度 约 1.2 m, 且 出 现 了 深度 为 0.7 m 左右 的 拉 
裂 损伤 区 。 由 图 9 可 知 ,从 监测 点 的 不 同位 置 分 析 ， 
在 隧道 据 进 的 瞬间 ,隧道 不 同 部 位 所 释放 的 弹性 应 
变 能 有 所 不 同 ,其 中 隧洞 两 侧 边 墙 所 储存 的 应 变 能 
略 高 于 拱 项 处 存储 的 应 变 能 ,在 左 侧 边 墙 处 最 大 ,最 
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大 值 为 2372 MJ/m 在 拱 顶 处 为 2115 MJ/m’。 


nr 弹性 应 变 能 密度 /(MJ，。m ) 
: T 2 380 


2 400[ 一 地 
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= 2000 1 904 
1 666 
= 1 600r- 1 428 
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EN 1200r 952 
ES d 714 
$i 476 


0. 
时 间 /ms 


监测 点 mu 
图 9 弹性 应 变 能 密度 时 空 分 布 (Hoek HAE) 


Fig.9 Space-time distribution of strain energy density 
( Hoek rockburst ) 

2) 基 于 Russenes 和 Turchaninov 判 据 分 析 

通过 数值 计算 , 获得 4 引水 隧洞 横 剖 面 的 
Russenes 和 Turchaninov 岩 爆 判 据 和 弹性 应 变 能 密 
度 时 空 分 布 ,如 图 10 ~ 图 13 所 示 。 由 图 10 ~ 图 13 
分 析 可 知 ,采用 Russenes 判 据 对 隧洞 进行 岩 爆 倾向 
性 分 析 时 ,隧洞 掌 子 面 的 左边 墙 与 左 拱 脚 附近 处 最 
大 指标 值 约 为 1.2, 指 标 值 大 于 0.55 的 深度 约 为 1.4 
m; 右 侧 边 墙 与 右 拱 脚 附近 处 最 大 指标 值 约 为 1.1, 指 
标 值 大 于 0.55 的 深度 约 为 0.9 m; 由 图 11 还 可 知 ， 
隧洞 左右 两 侧 边 墙 以 及 拱 顶 处 的 弹性 应 变 能 分 别 为 
2 420 MJ/m 、2 326 MJ/m 和 2 226 MJ/m 。 采 用 
Turchaninov 判 据 分 析 时 ,隧洞 掌 子 面 的 左 侧 边 墙 与 
左 拱 脚 附近 处 最 大 指标 值 约 为 1.4, 指标 值 大 于 
0.8 的 深度 约 为 2. 1 m; 右 侧 边 墙 与 右 拱 脚 附近 处 最 
大 指标 值 约 为 1.2 ,指标 值 大 于 0.8 的 深度 约 为 1.7 
m; 由 图 13 还 可 知 ,隧洞 左右 两 侧 边 墙 以 及 拱 顶 处 
的 弹性 应 变 能 分 别 为 2 334 MI/m' ,2 286 MJ/m 和 
2128 MJ/m’, 

通过 对 比 Hoek 判 据 .Russenes 判 据 和 Turchani- 
nov 判 据 对 岩 爆 倾向 性 的 分 析 结 果 可 知 ,三 者 数值 计 
算 结 果 均 表明 在 隧洞 关键 部 位 边 墙 处 , 左 侧 边 墙 相 
比 右 侧 边 墙 更 易 发 生 岩 爆 , 现场 实际 情况 则 表明 隧 
洞 掌 子 面 左 侧 边 墙 发 生 了 严重 岩 爆 , 而 右 侧 边 墙 并 
未 发 生 岩 爆 '2 ;但 是 ,在 预测 岩 爆发 生 的 位 置 方 面 ， 
Turchaninov 岩 爆 判 据 预 测 隧洞 掌 子 面 的 左 侧 边 墙 岩 
爆 位 置 更 靠近 隧洞 掌 子 面 的 墙角 部 位 , 因此 ， 
Turchaninov 判 据 与 Hoek 判 据 、Russenes 判 据 相 比 
较 , Turchaninov 判 据 预测 岩 爆 发 生 的 位 置 与 现场 情 
况 更 符合 ;在 预测 岩 爆 产生 的 破坏 深度 方面 ,三 者 数 
值 计 算 结 果 均 表明 对 隧洞 左 侧 边 墙 预 测 深度 较为 接 
近 ,Russenes 判 据 预 测 的 深度 略 小 于 Hoek 判 据 和 


Turchaninov 判 据 ,该 部 位 实际 岩 爆 发 生 的 破坏 深度 
约 为 2 m, 


ee SCC e CO C e t 
—i51muo-ooo-i 


10 Russenes 岩 爆 判 据 分 布 


Fig. 10 Russenes rockburst criterion distribution 
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SPig. 11 Space-time distribution of strain energy density 
< (Russenes rockburst ) 


图 12 Turchaninov 岩 爆 判 据 分 布 

Fig. 12 Turchaninov rockburst criterion distribution 

3) 基 于 ERR 和 LERR 指标 分 析 

由 于 目前 尚 不 明确 ERR 与 岩 爆 倾向 性 之 间 的 
关系 ,因此 无 法 依据 ERR 六 值 来 评估 岩 爆 的 倾向 性 
和 量化 岩 爆 的 烈度 等 级 。 但 ERR 方法 适用 于 施工 
方案 的 比 选 研 究 , 亦 可 总 体 上 描述 岩 爆 的 破坏 位 置 、 


范围 和 深度 。 例 如 ,运用 ERR 指标 ,采用 销 爆 法 对 
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隧洞 进行 全 断面 和 台阶 法 开 挖 ,所 得 隧洞 槛 剖面 内 
ERR 分 布 状 态 如 图 14 所 示 。 由 图 14 可 知 , 采 用 全 
断面 法 和 上 下 台阶 法 开 挖 后 计算 得 到 的 ERR BJf 4) 
别 为 2.54 x10” MJ/m? fll 2.34 x10 ”MJ/m ,这 表明 了 
采用 台阶 法 开 挖 隧道 诱发 岩 爆 的 风险 以 及 岩 爆 烈度 
明显 小 于 全 断面 法 开 挖 。 


HEr gn P 密度 /(MJ * m°) 
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图 13 弹性 应 变 能 密度 时 空 分 布 (Turchaninov £4) 


Fig.13 Space-time distribution of strain energy density 


( Turchaninov rockburst ) 
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(a) 台阶 法 开 挖 
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(b) 全 断面 开 挖 
图 14 ”隧洞 横 剖面 内 ERR 分 布 状态 (单位 :MJ/m’) 
Fig.14 ERR distribution of tunnel 


transverse section ( unit; MJ/m:’ ) 
通过 数值 计算 获得 的 隧洞 横 剖 面 内 LERR 和 弹 
性 应 变 能 密度 时 空 分 布 状态 ,如 图 15 ~ 图 16 所 示 。 
由 图 15 ~ 图 16 分 析 可 知 ,在 隧洞 掌 子 面 左 侧 边 墙 和 
右 侧 拱 腰 处 , 围 岩 局 部 能 量 释 放 最 多 , 且 分 布 范围 较 
广 ,LERR 最 大 值 为 2.91 x10? MJ/m 。 运 用 LERR 
判 据 指标 在 预测 隧洞 掌 子 面 左 右边 墙 及 拱 腰部 位 产 
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生 的 围 岩 能 量 释 放 率 时 ,隧洞 学 子 面 右 侧 边 墙 和 拱 
腰部 位 略 大 于 隧洞 掌 子 面 左 侧 边 墙 , 现场 实际 情况 
则 表明 隧洞 左 侧 边 墙 部 位 发 生 了 兰 爆 , 右 侧 边 墙 的 
混凝土 喷 层 出 现 了 膨胀 开裂 ( 见 图 17) ”。 
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15 ”隧洞 横 剖面 内 LERR 分 布 状 态 ( 单 位 :MJ/m ) 
Fig.15 LERR distribution of tunnel 
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- di HENEE uc. 
图 17 4* 引 水 隧洞 右 拱 肩 喷 层 鼓 胀 开裂 ”” 
和 
4) 基 于 围 岩 强度 比分 析 
由 图 18 ~ 图 19 分 析 可 知 : 隧 洞 最 大 岩 爆 坑 位 于 
左 侧 边 墙 及 拱 脚 处 ,这 与 现场 情况 接近 ,最 大 岩 爆 坑 
的 深度 约 为 5.6 m; 由 隧洞 现场 破坏 形状 可 知 ,采用 
围 岩 强度 比 作 为 岩 爆 判 据 的 数值 模拟 结果 与 实际 岩 
爆 坑 的 形状 ( 见 图 20) 基本 吻合 ,也 验证 了 围 岩 强度 
比 作为 岩 爆 判 据 的 合理 性 ,这 也 进一步 表明 洞 室 周 转 
岩 体 破裂 演化 的 平缓 加 速 历程 亦 是 围 岩 内 部 能 量 不 
断 积 聚 、 耗 散 的 过 程 ,应 力 值 的 突变 导致 围 岩 大 小 主 
应 力 差 变 大 ,使 得 能 量 积聚 增加 ,应 力 产生 高 度 集中 起 
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放 能 量 导 致 岩 体 出 现 大 变形 与 岩 爆破 坏 。 结 合 图 18 ~ 
图 19 可 知 ,隧洞 掌 子 面 右 侧 边 墙 及 拱 腰部 位 发 生 了 以 
剪 切 破坏 为 主 的 突 发 性 破坏 ,而 现场 实际 情况 中 该 位 
置 的 混凝土 喷 层 也 出 现 了 膨胀 开裂 ( 见 图 17 ) 。 
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图 18 围 岩 强度 比 


Fig. 18  Strength-stress ratio of surrounding rock 


图 19 岩 爆 模拟 示意 


Fig.19 Sketch of rockburst simulation 
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图 20 MAARRE 
Fig.20 Sketch of rockburst areas in situ 
5 ) 综合 分 析 

本 研究 在 数值 计算 的 基础 上 ,采用 Hoek 判 据 、 
Russenes 判 据 、Turchaninov 判 据 、ERR 指标 、LERR 
各 标 和 围 岩 强 度 比 6 个 指标 ,对 锦 屏 开 级 水 电站 4* 

引水 隧洞 的 岩 爆 倾向 性 进行 了 预测 评估 。 
Hoek 判 据 .Russenes 判 据 和 Turchaninov 判 握 三 
者 数值 计算 结果 均 表 明 , 引 水 隧洞 左右 两 侧 边 墙 发 
生 了 岩 爆 , 且 在 引水 隧洞 左 侧 边 墙 较 右 侧 边 墙 更 容 
易 发 生 岩 爆 , 现 场 实际 为 隧洞 左 侧 边 墙 出 现 了 严重 
岩 爆 , 右 侧 边 墙 并 未 发 生 岩 爆 ; 三 者 数值 计算 结果 表 
明 预 测 隧 洞 左 侧 边 墙 岩 爆破 坏 深度 为 1.5 mm~2.1m， 
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该 部 位 现场 实际 岩 爆 破坏 次 度 约 为 2m。 从 能 量 的 
角度 采用 ERR 指标 和 LERR 指标 进行 隧洞 岩 爆 倾 
向 性 分 析 可 知 ,隧洞 进行 全 断面 和 人 台阶 法 开 挖 时 ,人 台 
阶 法 发 生 知 爆 的 风险 以 及 岩 爆 烈度 明显 小 于 全 断面 
法 。 在 隧洞 堪 侧 边 墙 和 右 侧 拱 腰 位 置 处 围 岩 局 部 能 
基 释 放 区 域 最 多 , 且 分 布 范围 较 广 ,LERR 最 大 值 为 
2.91 x10° MJ/m  。 采 用 围 岩 强度 比 指标 分 析 时 , 隧 
洞 最 大 岩 爆 坑 位 于 掌 子 面 左 侧 边 墙 及 拱 脚 处 ,这 与 
现场 情况 接近 ,最 大 岩 爆 坑 的 深度 约 为 5.6 m; 隧洞 
掌 子 面 右 侧 边 墙 及 拱 腰 部 位 的 混凝土 喷 层 出 现 了 脱 
KFR. 


采用 上 述 6 个 判 据 指标 ,从 多 个 角度 预测 评 佑 
了 锦 屏 工 级 水 电站 4 引水 隧洞 的 兰 爆 倾向 性 (ERR 
判 据 指标 仅 进 行 了 开 挖 方案 的 比较 分 析 ) ,各 指标 的 


CN 基于 动力 数值 分 析 , 结合 Hoek 判 据 、Russenes 
判 据 .Turchaninov 判 据 、ERR 指标 .LERR 指标 和 围 
A838 8E EG 6 个 指标 ,并 以 锦 屏 开 级 水 电站 4 引水 隧 
洞 六 工程 背景 ,对 岩 爆 倾向 性 进行 了 定性 判断 ,得 到 
了 名 下 结论 。 

X) 在 预测 岩 爆 发 生 的 位 置 方 面 ,Turchaninov 岩 
爆 判 据 预 测 隧 洞 掌 子 面 的 左 侧 边 墙 岩 爆 位 置 更 靠近 
隧洞 掌 子 面 的 墙角 部 位 ,Turchaninov 判 据 预测 宕 爆 
发 生 的 位 置 与 现场 情况 更 符合 。 

2) 采 用 ERR 指标 对 开 挖 方案 进行 分 析 时 ,台阶 
法 开 挖 发 生涯 爆 的 风险 以 及 岩 爆 烈度 明显 小 于 全 断 
面 法 开 控 。 采 用 LERR 指标 进行 分 析 时 ,在 隧洞 左 侧 
边 墙 和 右 侧 拱 腰 位 置 处 围 岩 局 部 能 量 释放 区 域 最 多 ， 
且 分 布 范围 较 广 ,LERR 最 大 值 达到 2.91 x 10" MJ/m 。 

3) 采 用 围 岩 强度 比 指标 分 析 时 ,隧洞 掌 子 面 右 
侧 边 墙 及 拱 腰部 位 发 生 了 以 剪 切 破坏 为 主 的 突 发 性 
破坏 ,而 现场 实际 情况 中 该 位 置 的 混凝土 喷 层 也 出 
现 了 膨胀 开裂 。 
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